arXiv:1509.01635v3 [cond-mat.mes-hall] 16 Nov 2016 


Prepared for submission to JHEP 


Topology of Fermi Surfaces and Anomaly Inflows 


Alejandro Adem“ Omar Antolm Camarena" Gordon W. SemenofI* Daniel Sheinbaum'^ 

Department of Mathematics, University of British Columbia, 1984 Mathematics Road, Vancouver, 
Canada V6T 1Z2 

^Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, 
Vancouver, Canada VST IZl 

E-mail: adem@math.ubc.ca, oantolin@math.ubc. ca, gordonws@phas .ubc. ca, 
dshein@math.ubc.ca 

Abstract: We derive a rigorous classification of topologically stable Fermi surfaces of non¬ 
interacting, discrete translation-invariant systems from electronic band theory, adiabatic 
evolution and their topological interpretations. For systems on an infinite crystal it is 
shown that there can only be topologically unstable Fermi surfaces. For systems on a half¬ 
space and with a gapped bulk, our derivation naturally yields a A-theory classification. 
Given the d — 1-dimensional surface Brillouin zone of a d-dimensional half-space, our 
result implies that different classes of globally stable Fermi surfaces belong in K~^ (Xg) for 
systems with only discrete translation-invariance. This result has a chiral anomaly inflow 
interpretation, as it reduces to the spectral flow for d = 2. Through equivariant homotopy 
methods we extend these results for symmetry classes AI, AH, C and D and discuss their 
corresponding anomaly inflow interpretation. 
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1 Introduction 

Classification of locally stable Fermi surfaces and topological phases of matter using K- 
theory was first introduced in the pioneering work of Hofava [1] and Kitaev [2]. The present 
work generalizes to the infinite dimensional Hilbert space setting and global character the 
analysis of Hofava, though, through a rather different route. We do not restrict ourselves 
to analyzing the local character of perturbations around a chosen fixed-point of the Fermi 
surface. However, for the sake of completeness, we illustrate Hofava’s analysis in this 
setting. The main novelty in our approach is the fact that, to our knowledge, we are the 
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first to obtain directly a if-theoretic classification for global stability of Fermi surfaces 
without resorting to physically-unjustified mathematical procedures, as it has been done 
previously in the literature. The chiral anomaly connects to topology via the Atiyah- 
Patodi-Singer index theorem [3] and the notion of spectral flow. This topological invariant 
comes out in our construction. Edge states of the Integer quantum Hall effect (IQHE) were 
already given a spectral flow interpretation by Fukui and collaborators [4]. Here we put 
their result in the wider framework of A-theory and extend to other symmetry classes. 

Ryu et al. [5],[6] and Witten [7] among others, have developed a formalism for studying 
gapped bulk systems through classifying their gapless edge states and the bulk-boundary 
correspondence, where the edge states have an anomaly inflow interpretation [8], [9]. This 
paradigm is physically similar to ours, however these arguments rely upon taking a low- 
energy limit of a lattice Hamiltonian to obtain an effective relativistic field theory [10] such 
as Chern-Simons for the IQHE, with a correspondence between topological field theories 
and symmetry protected topological phases [7], [11]. The present work does not rely on 
such approximations, but we do have other strong assumptions such as discrete translation- 
invariance. 

The results presented below show that systems on an infinite crystal cannot posses 
globally stable Fermi surfaces, up to homotopy. The proof presented here is, to our knowl¬ 
edge, new and does not rely on the usual arguments of fermion doubling [12]. It is also quite 
simple from a mathematical point of view. We proceed to study systems on a half-space 
X [0, oo). As discrete translation-invariance is broken by the surface x {0} only 
a reduced version of Bloch’s theorem holds. Nevertheless this opens up the possibility of 
having surface states and so called bulk states become part of the continuous spectrum of 
our system. This, together with a gapped bulk condition will play the central role in the 
derivation of the A-theoretic classification of globally stable Fermi surfaces for systems on 
a half-space. 

Issues of aperiodicity, disorder and unitary-antiunitary ambiguities in representations 
of symmetry operators are considered in a more general framework in [13], [14] and will not 
be dealt with in this work. However, neither of [13], [14] attempt to classify Fermi surfaces. 
We work in the so-called first quantization formalism where we construct the ground state 
of our systems employing the single-particle picture and filling the lowest single-particle 
states using Pauli’s exclusion principle. 

2 Momentum-space classification with only discrete translation-symmetry 

We proceed to derive a classification for both globally topologically stable Fermi surfaces 
for systems on an inhnite crystal in subsection 2.1 and half-spaces x [ 0 , 00 ) in sub¬ 

section 2.4. To put our classification in the adequate context, we include in subection 
2.3 a simplified version of Hofava’s classification for locally stable Fermi surfaces in the 
single-particle setting, as well as its limitations for modeling globally stable surfaces. By a 
globally topologically stable Fermi surface we mean that its associated gapless system can¬ 
not become gapped (have an empty Fermi surface) as long as certain conditions, carefully 
detailed below, are satisfied. Otherwise the Fermi surface is referred to as topologically 
globally unstable. 

All systems considered here are sometimes treated in the literature in the second quan¬ 
tization formalism, constructed in a fermionic Fock space of particle creation/annihilation 
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operators. As stated previously, we work under the single-particle picture and we have 
adapted Hofava’s local classification into first quantization form, which allows us to en¬ 
compass the different systems with the same formalism and also permits a simpler inter¬ 
pretation of the operator theoretic aspects of the problem, as in [13], [14], 

2.1 Infinite crystal 

Let X be the Brillouin zone of a non-interacting fermionic system on an infinite crystal, i.e. 
W^. For each crystal momentum ^ G X there is a Bloch Hamiltonian operator acting 
on a complex separable Hilbert space J^. %{k) is self-adjoint and there exists a complete 
orthonormal basis of such that 


%{%) |4>n(^)) = £n(k) |4>n(fc)) , 

(2.1) 

n{k) = (zV - + H(r), 

(2.2) 

V{f^ = V{r + R) 

(2.3) 


where the A:-dependent eigenvalues en{k) are known as energy bands and ii is a lattice 
vector. 

The Fermi surface [15] of a non-interacting fermionic system corresponds to the set of 
points in the Brillouin zone for which one or more energy bands are equal to the Fermi 
energy £p 


G X I {T-L{k) — £pl) \^n{k)) = 0 , for some n| (2-4) 

Hence the operator of interest is i'H{k) — £pl) and we wish to study its kernel, which has a 
physical interpretation. We have the freedom to choose our energy scale such that = 0 
and we only write %{k) instead. If we restrict ourselves to the n-th energy band £nik), 
the corresponding surface embedded in the Brillouin zone is called the n-th branch of the 
Fermi surface, the Fermi surface being the union of all such branches. We assume that 
none of our Bloch Hamiltonians are infinitely degenerate at the Fermi energy, thus 

dim Ker%(k) < oo 'ik€X. (2.5) 

The work of Atiyah and Singer [16] is central to our classification. To enable the use of the 
machinery developed in [16] we must impose that our Hamiltonian operators be bounded. 
This is a standard trick employed when analyzing topological properties of Hamiltonian 
and Dirac operators, which is usually done by employing the topology induced by the Riesz 
norm, defined by imposing that the transformation 

n{k)e^n{k){i+ n{k)^'Hik ))--2 ( 2 . 6 ) 

be an isometry^ making unbounded operators bounded [17]. Thus, we could avoid the 
unnatural bounded condition by working in the so called Riesz topology, but for the sake 
of simplicity we shall stick with the operator norm topology for the time being. We give 
the details of the equivalence of both treatments in Appendix A. We list the assumed 
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properties of our Bloch Hamiltonian in mathematical terms: 

• T-LQt) is bounded, 

• dim Ker'H{k) < oo , 

• T-LQt) is self-adjoint. 

Operators satisfying these conditions are known in the mathematical literature as bounded 
self-adjoint Fredholm operators [16], where the set of bounded Fredholm operators [18] 
are bounded operators in with finite kernel and cokernel. We denote the space of all 
complex bounded Fredholm operators and bounded self-adjoint Fredholm operators 

as 

2.1.1 Path-components of the space of bounded self-adjoint Fredholm opera¬ 
tors 

Let us endow F{J^) with the norm topology, the most natural topology for bounded 
Fredholm operators (see Appendix A). With this choice of topology has 3 disjoint 

path-components, operators which are essentially positive denoted essentially 

negative and the remaining path-component Care must be given to the 

subtleties of what it means to be essentially positive/negative. A self-adjoint operator B is 
positive if 

(2.7) 

and essentially positive if 3 iF c J'F' with codim F < cq such that B is positive on F and 
BF = F [16], [19]. The definition of essentially negative is analogous. Depending on which 
path-component our Hamiltonian belongs to we will obtain very different results. 

2.2 Adiabatic evolution of systems on an infinite crystal 

77(fc) is a continuous function of k, so its image lies in one of the 3 disjoint path-components 
of With this in mind, we wish to put in the same equivalence relation all systems 

which can be adiabatically evolved [20] into one another. As usually done in the literature, 
we employ a homotopy model of adiabatic evolution, i.e. V 77o, 77i : X —)■ we have: 

77o ~ 77i 3g :Xx I ^ continuous, 

g{k,0)=no{k), g(fc,l)=77i(fc) VfcEX. (2.8) 

There is a rigorous notion of adiabatic evolution in condensed matter systems given by 
Avron et al. [20] with applications to the IQHE, and it has been generalized to systems 
without a gap condition [21]. We emphasize that eq.(2.8) is a sufficiently good approxima¬ 
tion to [20] only when we restrict the space of allowed Hamiltonians and is in general not 
true for arbitrary ones. We also remark that for metals the notion of adiabatic evolution 
should be understood in the context of Landau’s Fermi liquid theory [22] as in [1],[23],[24]. 

Let us denote by (T{T-L{h)) the spectrum of T-L{k). On an infinite crystal Bloch’s theorem 
induces the following periodicity condition 

4>nik,r-\-R) = 4>nik,r), (2.9) 

were again R is an element of the direct lattice that represents the crystal. Let us remember 
that the labeling of solutions by the crystal momentum k is due to the unitarity of the 
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translation operator and Bloch’s theorem [25] and since k k has real components. 
Furthermore, if l^(r) = V{f + R), is also bounded, then ^{{k) has compact resolvent, 
implying that the continuous spectrum of T-L{k), adT-Lik)) is empty. Thus the spectrum 
of 7i{k) coincides with its eigenvalues. Since T-l^k) is bounded from below, it has a finite 
number of valence bands (negative eigenvalues for each k) [25], i.e. R^k) G 
after applying a Riesz transformation (Appendix A). Thus, all types of Fermi surfaces for 
systems on an inhnite crystal, separated by equivalence relation (2.8) are given by the set 
of homotopy classes of maps [X, J^^“(J^)]. However, Atiyah and Singer proved in [16] that 
is homotopically equivalent to a point via the simple homotopy 

'Ht{k) = {l-t)U{k) + mo{k), (2.10) 

where T-Lo{k) is any gapped bounded Bloch Hamiltonian V A: G X. This implies 

[X,Rfd^)] = *, ( 2 . 11 ) 

* denoting a single, i.e. trivial class. Hence, all systems as described above can be con¬ 
nected to a gapped system through adiabatic evolution. This leads us to the following 
conclusion: all systems in symmetry classes A, AI, and AH, on an infinite crystal, only 
have globally topologically unstable Fermi surfaces, up to adiabatic evolution. Class A refers 
to systems whose only symmetry is discrete translation-invariance in the Altland-Zirnbauer 
classihcation [26]. The validity for the remaining symmetry classes AI and AH will be 
discussed in section 4. 

This result was expected by arguments similar to those of fermion doubling on lattice 
gauge theories [12], by the Nielsen-Ninomiya theorem, where the Hamiltonians under con¬ 
sideration are matrices. For a recent discussion of the role of the Nielsen-Ninomiya theorem 
in condensed matter see [27]. Note that though our result yields sufficient conditions to 
guarantee the topological instability of the Fermi surface of a system on an infinite crystal, 
it does not have anything to say about necessary conditions. In the argument presented 
above we used the fact that for the infinite crystal, the spectrum of T-l{k) is discrete and 
bounded from below, which guarantees that ^{{k) G These are not necessary 

conditions as there can also be unbounded Dirac Hamiltonians, which do not belong to 
but still correspond to a trivial class and have a globally topologically unstable 
Fermi surface. We also remark that the topology of the Brillouin zone was irrelevant for 
this result since what produces the triviality is that X'™(=^) is contractible [28]. 

Finally, even though the physical result was expected, the method employed here is 
new to our knowledge, and perhaps the simplest and most straightforward method to prove 
such a result. It does not depend on the dimension of the system and does not depend on 
taking a low-energy limit of a lattice Hamiltonian [10]. 

2.3 Local stability of Fermi surfaces 

Hofava [1] originally developed a classification of locally stable Fermi surfaces, that is, 
pieces of the Fermi surfaces which are robust to small perturbations near a given ko on the 
Fermi surface. Later on Matsuura et al. [23] and Chiu et al. [24] extended Hofava’s work to 
include more symmetry classes. Both [23],[24] clarify that they restrict to a neighborhood 
around a point because of the similar arguments to the above ones concerning fermion 
doubling. From our construction it would seem that restricting to a subspace of X does 
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not yield a different result since (2.11) is true for any choice of space X. What we must do 
is also restrict the set of allowed homotopies. Let us denote one of the path-components 
of the Fermi surfaces by S and let dim Tj = d — p — 1. One looks at the p + 1-transverse 
directions to S in X around a chosen point /cq G S. Hofava’s resolution, for which we 
present a simplified version, adapted to the infinite dimensional formalism, is to consider 
the boundary of a ball dBrg{ko) C X, which does not intersect S. The boundary is a 
topological sphere S^, whose radius ro = ro{ko) depends on the choice of fco G S. As long 
as we impose that the adiabatic evolution of the system does not generate a zero eigenvalue 
for any T-L{k), with A: G S^, we have restricted homotopy classes of maps 

[SP,TV{Jf)]^ Vectc{^P). ( 2 . 12 ) 

where TV{^) C denotes the space of invertible Hamiltonians [13] with empty 

continuous spectrum, also called gapped, i.e. 

dimKer%{k) = 0 'i'H{k) G TV{J^) , (2.13) 

ac{n{k)) = <j). (2.14) 


These operators were widely studied in [13], [14], [29],[30]. Fectc(S^) is set of isomor¬ 
phism classes of complex vector bundles with base space [31]. Equation (2.12) can 
be derived by repeating a simplified version of the analysis in [13], section 10. The ar¬ 
guments presented there follow through, except that in our case, we must restrict the 
parameter space for our families of Bloch Hamiltonians to dB^^iko) ~ instead of the 
whole Brillouin zone X. If we are willing to lose information on the homotopy classes but 
gain computability via a stabilization (see [13], [31] for details, or [29],[30] for an approach 
without stabilization) we obtain the reduced A-theory group 


A(SP) 


Z p even 
0 otherwise. 


(2.15) 


Here K{SP) denotes the reduced Grothendieck completion of Vectc{S^), the semi-group 
of isomorphim classes of complex vector bundles with base space [31]. 


2.3.1 Limitations of local stability 

In order to avoid our hrst conclusion about the topological triviality of Fermi surfaces on an 
infinite crystal, we had to restrict the set of allowed homotopies in the analysis presented 
above. There is no natural nor universal procedure to make such a restriction for all /c G X 
as there is no analogous quantity that would play the role of the gap in the bulk spectrum 
for topological phases for the whole Brillouin zone. Thus, we would have to look at the 
particularities of each chosen point /cq G S to define an adequate S^. In particular this 
construction is not well suited for modeling globally stable gapless edge modes observed 
at the boundaries of topological phases of matter. That being said, because of the role 
of TV{Jif) in Hofava’s analysis, under the right codimension d — p — 1, the local analysis 
emulates part of the global analysis of gapped topological phases (see [23], [24], [32] for a 
discussions and connections) but it strongly fails to account for all the involution fixed- 
points induced by symmetries [24], which give rise to distinct topological phases [13]. 
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2.4 Systems on a half-space 

Let us now consider systems on the half-space x [0,oo), i.e. systems with a single 

boundary. Typically the boundary is a domain wall between two different systems (or 
vacuums). Let us denote the directions parallel to the the boundary by fj|. Since we no 
longer have the periodicity of the potential in the direction perpendicular to the boundary 
(denoted as r_L), we can only apply a partial Bloch theorem, in the directions parallel to 
the boundary. We set the boundary at r_L = 0. Our wave functions obey 

i^a(k\\,f\\,rA_) = e*^llOl(/)„(fe||,f^|,r_L), (2.16) 

^R|l = l</>a(fc||,r_L)) . (2.17) 

Where i?|| is a surface lattice vector and is the associated discrete translation 
operator. Let us remark that, similarly to the infinite crystal case, the labeling of solutions 
by fey comes from the unitarity of and hence we can restrict fey to be in X^, where 
Xs is the Brillouin zone associated to the surface reciprocal lattice. There are different 
boundary conditions that can be imposed at r± = 0, depending on the system under 
consideration. In general one must match a solution with the same functional form as an 
infinite crystal solution (i.e with no boundary) to solutions of the free Schrodinger equation 
(or more generally to the solutions of a Schrodinger equation in a different crystal), which 
decays exponentially as r_i_ becomes positive and we leave the crystal [15], [33], [34]. The 
infinite crystal-like piece of '(/’□.(fey) has an factor, where k± originally represented 

the component of the crystal momentum corresponding to r±. If we had an infinite crystal, 
fe_L would be forced by unitarity of to be real. Nonetheless, since the boundary breaks 
the discrete translational symmetry in the r_L-direction, k± is no longer a good quantum 
number and it can take complex values as well (unlike the components of fey, which must 
be real by construction) [15], [33], [34]. These complex values of fe_L allow for the existence of 
new kinds of solutions described below. We thus transform the single-particle Schrodinger 


equation into one for each fey of the form 


77(fcy) \(l)a{k\\,r±)) = ea(fcy) \4>a{k\\,r±)) , 

(2.18) 

77(fey) = + (fVy - fey)2 + V{f\\,r^) , 

(2.19) 

^(pi ) ^±) = ^ (pi + -R||) ^±) • 

(2.20) 


The index a = a{k±,m) behaves similarly to the principal quantum number n of the hy¬ 
drogen atom, i.e. ea{k\\) can either be part of the continuous spectrum of ?7(fey), (Tc(?7(fey)) 
or part of the discrete spectrum 0-^(7^ (fey)). 

There are 3 types of solutions to the associated PDE problem, represented in figure 1. 
Solutions with real fe_L {Im{k±_) = 0) are called bulk states and correspond to truncating at 
the boundary the periodic solutions of the original infinite crystal. These can have any real 
value of fe_L and thus, they are part of ac{'H{k^\)). Note that though the PDE problem has 
a solution, the eigenvalue equation has no solution in for eo(fey) G (Tc(77(fey)), similar 
to the scattering states of the hydrogen atom [35]. Solutions with Im{k±) / 0 not only 
decay exponentially outside the crystal but also have an exponentially decaying envelope 
in the r_i_-direction, towards the crystal 

4>a{k\\,r\\,r±) « e"l^”'(^-L)^-Llua(fey,f||) (2.21) 
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(I) 



( 11 ) 



(III) 



Figure 1: Real part of the wave function of the distinct types of states on a half-space as 
a function of r±. All solutions decay exponentially outside the crystal. II is a bulk state 
which is periodic inside the crystal. Ill Surface state, decays very rapidly away from the 
boundary, fill Resonance state, the decay eventually stops and the solution is again 
periodic. Dashed lines represent the envelope of the solutions. 


where Ua{k\\,f\\) is a bounded and periodic function of r^|. These solutions are called 
surface states [15],[33], [34]. Matching such solutions to vacuum solutions at the boundary 
forces k± to come in discrete values, implying that surface states correspond to £aik\\) G 
ad{'H{k^^)) [15],[34]. m G Z is the number that arises from having Bloch’s theorem in d — 1 
dimensions. Note that surface states do belong in as they are normalizable. Finally, 
there can be solutions which join a surface state with a bulk state, called surface resonance 
states and also form part of ac.{fl-L{k\\)). This joining is a limiting process in which a 
bulk/surface state becomes a surface/bulk state at some A:|| G X^, i.e. Ua{k\\,r\\) continues 
to be periodic but stops being bounded by 

2.5 Adiabatic evolution of systems on half-space 

If we impose that the bulk be gapped, the Fermi energy Sp must he between bulk states, 
i.e. Ep ^ crc('R(A:||)). Setting Sp = 0, this implies that ac{'H{k\\)) n / (j), corresponding 
to the bulk valence band nearest to the Fermi energy and so 'R(A:||) ^ Similarly 

ac{'H{k^\)) n M"*" / corresponding to the bulk conduction band nearest to the Fermi 
energy, which implies that H(fe||) ^ Thus, H(fe||) G as long as the 

bulk remains gapped. Globally stable Fermi surfaces of systems on a half-space are thus 
classihed by the set of homotopy classes [X^, T]f“(^)]. 

At this point let us remark that adiabatic evolution for systems with metallic surfaces 
refers to a bulk adiabatic evolution, i.e. bulk states below Ep cannot tunnel directly to bulk 
states above the Fermi energy under adiabatic evolution, as defined in eq.(2.8). Nonetheless 
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bulk/surfaces states are allowed to tunnel to surface/bulk states through surface resonances. 

The set of skew-adjoint Fredholm operators J-{^) = is trivially homeomor- 

phic to {^) in the norm topology, and so are T"™ (=^) and T"* (J^) , the corresponding 
non-trivial component of Atiyah and Singer proved in [16] that ~ 

the loop space of the set of Fredholm operators Combining this with Atiyah’s proof 

[18] that is a classifying space for complex A-theory and the suspension isomor¬ 

phism [28] we get 

[x„JT(^)] « 

^ [X„11J^(^)], 

« [SXs,T{je)],, 

-K-\Xs), ( 2 . 22 ) 

where K~^{Xs) ss K{SXs^^) denotes the reduced A-group and SXg^^ is the suspension of 
our surface Brillouin zone with a disjoint point attached [31]. Notice that Grothendieck’s 
completion popped out naturally from our construction! 

3 A-theory or vector bundles? 

Geometric constructions of fermionic topological phases on an infinite crystal with no 
boundary often lead to isomorphism classes of n-dimensional vector bundles over the d- 
dimensional Brillouin zone X, Vect'^iX). In general, passing from Vecf^fX) to A-theory 
involves a series of steps. First we must allow addition of isomorphism classes [Ai] © 
[A 2 ] = [Ai © A 2 ] and adding 1 -dimensional vector bundles (i.e. single bands), which forces 
us to consider isomorphism classes of vector bundles of any dimension over X, arriving 
at the semigroup VectciX). Then, we can either impose the stabilization equivalence 
relation [31], as we did in subsection 2.3 to obtain the reduced A(X) group or we can 
define another equivalence relation, defined on pairs ([Ai], [A 2 ]), ([A(], [A 2 ]) of isomorphism 
classes, turning VectciX) into the abelian group A(X) [31]. The latter process is called 
a Grothendieck completion or Grothendieck construction and is usually ignored in the 
physics literature, where most works compute the unreduced A(X) group directly to classify 
fermionic gapped phases, without justifying the transition to A-theory. As adequately put 
by Thiang [14], this has produced a conflation between the full A(X) groups and their 
reduced versions A(X), leading to discrepancies among some tables of gapped topological 
phases [14]. 

The mathematical literature on the other hand addresses the issue in different ways. 
As examples we can consider the work of De Nittis and Gomi [29],[30], which uses n- 
dimensional vector bundles and their characteristic classes to classify fermionic phases, 
without appealing to A-theory. In a different vein Freed and Moore [13] employ a process 
equivalent to a Grothendieck completion but emphasize that they are not aware of an 
apriori good physical motivation for its use. There are also examples in the C^-algebraic 
approach of Prodan and Schulz-Baldes [36] together with references therein. Prodan and 
Schulz-Baldes motivate the use of A-theory for lattice Hamiltonians in the so called tight- 
binding approximation [15], through its applications in the bulk-boundary correspondence 
when combined with cyclic cohomology. Though a classification using A-theory is not a 
direct logical consequence of their formalism, its use in their approach is certainly well 


- 9 - 


motivated and it employs ilT-theory of C^-algebras for both surface states and bulk states, 
with similar results as ours for the surface states in symmetry class A, employing the 
functor K~^ but on a real space algebra of surface observables. 

To our knowledge, our construction for systems on a half-space is the only one for which 
/f-theory is not put in by hand but rather arises naturally from the simple mathematical 
interpretations of adiabatic evolution and a gapped bulk spectrum for Bloch Hamiltonians 
on a half-space. We also do not restrict to Hamiltonians in the tight-binding approximation 
nor lattice systems, nevertheless our construction has the rather strong assumption of 
discrete translation-invariance on the boundary and the more common assumption of a 
gapped bulk spectrum, restrictions which are better handled in [36] . 

4 Adding symmetries 

Let us mention from the offset that we will only study systems on a half-space for which 
the boundary conditions preserve the symmetry class of the corresponding infinite crystal 
(bulk) Hamiltonian. Thus, both the continuous and discrete spectrum of our Hamilto¬ 
nians ac{'H{k^\)) respectively, inherit the same structure when a symmetry is 

implemented. 

4.1 Real and quaternionic structure on a complex vector space 

We must first digress and introduce the following mathematical definitions taken from [37]. 
A real structure on a complex vector space V is an anti-linear operator K, such that 
= I. In this case V ~ W (8 >m C, where W is vector space over M. A quaternionic 
structure on a complex vector space V is an anti-linear operator K, such that = —I. 
I, K, i and J = iK satisfy the quaternion relations. As will be shown below these structures 
will be induced on our Hilbert space by the corresponding symmetry operators and will 
modify the complex A-theory that appears in our classification into a different (but similar) 
extraordinary cohomology theory [28], that we will use to compute the distinct classes of 
globally topologically stable Fermi surfaces in a given symmetry class. 

4.2 Particle-hole symmetry 

The first symmetry we consider is the so called particle-hole symmetry. Systems with 
only particle-hole symmetry and discrete translation-invariance are denoted as symmetry 
classes C and D in the Altland-Zirnbauer classification [26] and the operator representing 
said symmetry is traditionally denoted by H. There are different choices of implementations 
of H one can make, such as whether it is a unitary or antiunitary operator [13],[14] and each 
one of them represents systems in condensed matter with different physical properties. Our 
choice of implementation is the one most commonly used throughout the physics literature 
[2], [23] [38], where our symmetry operator H satisfies 

HA(fc||) = -A(-fc||)H, (4.1) 

“i = -iB . (4.2) 

This implies that for each spectrum index a of 1-L{k\\), there exists a spectrum index /? of 
'A(—/c||), such that 

^\(l>a{k\\,r±)) = \(j)/3{-k\\,r±)) , (4.3) 

e„(fe||) = -e/3 (-fe||). (4.4) 
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In general /3 / a to respect the symmetry, however if £a{k\\) £ ^diT~L{k^\)) passes through 
an involution fixed-point /cji = —fey, then /3 may equal a. This is a subtle but non-trivial 
point as it permits the existence of some topologically non-trivial Fermi surfaces for all 
symmetry classes considered here (see subfigures 3II, 3III). Note that for particle-hole 
symmetric systems there are no J^^“-components due to equation (4.4). 

4.2.1 Class D 

If our system belongs to Class D, we must further assume that 

h2 = /. (4.5) 

Thus, H represents a real structure on our complex Hilbert space as defined above, 
which we denote as It also induces an involution t : Xg ^ Xg, r(fc||) = —k\\, making 
Xg a real space in the sense of [39]. The general strategy for extending the results of [16] 
is to employ equivariant homotopy [40], where the above implementation of the symmetry 
induces the action ?^(/c||) i—)• —H?^(—fc||)H. We wish to use a slightly different action, defined 
in [41], for which we must pass to the skew-adjoint operator i?^(/c||), with the action 

f?^(fc||) Ei'H{—^)E. (4.6) 

This subtle difference allows us to employ a result of Matumoto [41] 

[Xg,T{J^R)]c,-KR{Xg). (4.7) 

where KR is a different kind of if-theory, originally developed in [39] and C 2 denotes the 
cyclic group of order 2, with the associated action (4.6). Using Whitehead’s equivalence in 
conjunction with (4.7), it suffices to have plain homotopy equivalence on the fixed-points 
of the subgroups of C 2 to have the full equivariant equivalence, yielding 

\Xg,T{M’R)\c^ ^ KR-\Xg ). (4.8) 

The equivalence on the corresponding fixed-points is given by the Atiyah-Singer map em¬ 
ployed in [16], said map being equivariant under the chosen action. This group represents 
the topologically stable isomorphism classes of Fermi surfaces in symmetry class D. Af¬ 
ter seeing the implications particle-hole symmetry has for a{l-L{k\\)) it becomes clear that 
such a symmetry cannot be fully realized in condensed matter systems, nevertheless it is a 
powerful approximation to employ in the analysis of superconductors [42] where it is true 
for states near the superconducting gap A. 

4.2.2 Class C 

Now we move on to the case where 

“2 = -/ (4.9) 

also known as symmetry class C. For this case, the particle-hole symmetry operator in¬ 
duces a quaternionic structure on iE being the fourth generator. Let us denote this 
quaternionic structure by Jifg. Again we employ the same strategy as above, passing to 
skew-adjoint operators where, following [41] we denote by D 2 the cyclic group of order 2, 
which acts on by the different action 

i?7(fe||) ^-S'H(-fe||)H. (4.10) 
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Matumoto also proved the equivalence 


[Xs,TiJ^Q)]D,-KH{Xs). (4.11) 

KH is known as symplectic KR-th.eoiy, and we will follow the notation employed in [43]. 
Repeating the fixed-point argument as above, we conclude that stable isomorphism classes 
of Fermi surfaces in class C are given by 

[Xs,fiM’Q)]D, ~ KH-\Xs ). (4.12) 

We can compute KH~^{X) for any n entirely in terms of KR-gioups, since we have the 
isomorphism 

KH-^{X) « KR-^-^{X). (4.13) 

4.3 Time-reversal symmetry 

We shall implement time-reversal symmetry through an operator 0 [2], [38] 

0H(fe||) = H(-fc||)0, (4.14) 

Qi = -iQ. (4.15) 

Similarly to particle-hole symmetry, this implies that for each spectrum index a of H(fe||), 
there exists a spectrum index j3 of such that 

Q\Mk\\,r±)) = \4>0{-kii,r±)) , (4.16) 

Saikii) = epi-kii). (4.17) 

There are two classes in the Altland-Zirnbauer classification, that correspond to sys¬ 
tems with only discrete translation-invariance and time-reversal symmetry, namely classes 
AI and AH, whose properties we describe below. 

4.3.1 Class AI 

Systems in class AI must obey 

02 = /. (4.18) 

We now have the appropriate action directly on self-adjoint operators since 

n{k\i) ^ en{-k\\)e. ( 4 . 19 ) 

Looking again at the fixed-points of the subgroups of C 2 using the homotopy equivalence 

J-™(^) ~ (4.20) 

proven in [16], which holds for the fixed-points of all subgroups of C 2 and where is a 
=t=-represantation of the Cq Clifford algebra or its complexification, depending on the choice 
of subgroup [16]. This equivalence is equivariant with respect to (4.6), hence, together with 
(4.7) this implies that 

[X., JT(=4^h)]c2 ~ KR-'^{Xs) , (4.21) 

corresponds to isomorphism classes of stable Fermi surfaces in symmetry class AI. 
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4.3.2 Class All 

Finally, for stable Fermi surfaces in symmetry class AH, we must assume 

©2 = -/. (4.22) 

By repeating the same process for self-adjoint operators as we did for skew-adjoint operators 
in class C, since 

H(fc||) ^-0?^(-fc||)0 (4.23) 

and we obtain 

[Xs,J^ri^Q)]D, ~ KH-\Xs ). (4.24) 

We remark that both H and 0 are either real or quaternionic structures on our physical 
Hilbert space and that the difference between them lies in the set of fixed-points under (4.6), 
(4.10), (4.19), (4.23) in the space of Hamiltonian operators FI°'{Mq) or their 

action-induced skew-adjoint counterparts. If we repeat the above process for classes AI and 
AH but instead look at equivariant homotopy classes of the J^™-component, corresponding 
to systems on an infinite crystal with no boundary, we see that they are equivariantly 
contractible, proving our first stated conclusion in subsection 2.1. 

4.4 Relation to the work of Preed-Moore 

A much more general method was developed by Freed and Moore [13], who undertook the 
arduous task of laying out all possible representations of symmetries and how they modify 
the classification for gapped systems (for a C^-algebraic version see [14]). Unfortunately 
this method doesn’t admit a trivial extension to our construction. In a few words, this 
happens because for the general case, symmetries need the compact-open topology but the 
space of Fredholm operators is contractible in this topology. The pattern underlying these 
examples and the work in [13],[14] suggests that to include symmetries it is necessary to 
have a generalized twisted equivariant A-theory [44],[45], that reduces to the variants we 
have mentioned here. 

5 Analysis 

We shall compute different classes of topologically stable Fermi surfaces for the five cases 
discussed above (discrete translation-invariance, particle-hole symmetry with = ±7 and 
time-reversal symmetry with 0^ = ±7 ) for d-dimensional systems on a half-space, i.e. 
when the surface Brillouin zone ad— 1-dimensional torus or its involutive 

version for d = 1,2,3. We only present with some detail the computations of the 

different classes of Fermi surfaces for the only non-trivial case when d = 1-dimensional 
systems on a half-line and d = 2-dimensional systems on a half-plane (Fermi points). 

5.1 d = 1 systems 

The only non-trivial example for d = 1, = /cq, a point, is in symmetry class D, where 

KR-\ko) « KR{S^), 

^ KO{S^) 

« Z2 (5.1) 
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(I) (11) 
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^0 

0 -|- 0 -f- eo ^ 0 

-^0 
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Figure 2: d = 1 unstable and stable modes. 21 has no zero modes, where as 2II 
represents a topologically stable Majorana zero mode as an edge state, /cq is the 
0-dimensional surface Brillouin zone. 


The non-trivial element of Z 2 represents the so called Majorana chain developed by Kitaev 
[46], shown in subfigure 2II. The single-particle picture corresponds to a Bogoliubov quasi¬ 
particle and we have taken the mean field approximation to electron-electron and electron- 
phonon interactions [42]. Note that the validity of adiabatic evolution (Fermi liquid theory) 
for d = 1-systems is strongly limited. 

5.2 d = 2 systems 

For d = 2, our Brillouin zone Xg = is a circle. When we only have discrete translation- 
invariance. 


K-\S^) ^ k{S^) 

«Z. (5.2) 

The topological invariant that determines the corresponding class of a Fermi surface in 
K~^(S) ~ is known as the spectral flow [3]. This non-trivial state is shown 

in subfigure 31. The spectral flow had been employed previously to study chiral anomalies 
[47]. For this particular example, the physical interpretation of the spectral flow of our 2- 
dimensional system corresponds to either bulk valence/conduction bands (positive/negative 
flow) having surface resonances, which become surface states, eventually cross the Fermi 
energy and become again surface resonance states in the bulk conduction/valence bands 
respectively. This amounts to a charge C/(l)-chiral anomaly inflow [7], [27] when an electric 
field is applied to the system. This is because an electron in a bulk valence/conduction 
band state becomes localized at the surface before eventually becoming a bulk conduc¬ 
tion/valence state and it would appear that the boundary has gained/lost electric charge, 
when it is actually incoming/leaking from/into the bulk. 

Realizations of such systems are surface states of the IQHE [38], once magnetic trans¬ 
lations have been taking into account. Fukui and collaborators had already given a spectral 
flow interpretation to the IQHE [4], and further relate the spectral flow to a bulk Chern 
number but did not relate it to ET-theory as we have done in this work. 

If we now consider a 2-dimensional system with = /, the involution r(A:||) = — A:|| is 
equivalent to having the Brillouin circle embedded in C, with r being complex conjugation. 
It has both {—1,1} as its fixed-points and following [43] we denote our Brillouin circle 
with involution as where is the {p + q — l)-dimensional sphere embedded in 
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Figure 3: d = 2 Topologically stable Fermi points. Solid regions and dashed lines belong 
to the continuous spectrum ac{'H{k\\)) (bulk states and resonance states respectively). 
Solid lines to (surface state). + (—) signs label the sign of the slope at 

crossings. 31 corresponds to class A (IQHE) with a net spectral flow = 1. 3II corresponds 
to a spectral flow of 2, that may occur for both classes C and D at non hxed-points in 
3III represents two Majorana zero modes in class D. 3IV represents a helical Dirac 
fermion edge mode (SQHE). 


= RP 0 iM*?, Thus we obtain 


KR-\S^’^) ^ KO-^®KO 

RiZeZa, (5.3) 

where KO^ denotes the real (trivial involution) j-th iT-group of a point. For more details 
on how to compute KR^{X.) we refer to [43], section 3. The Z-factor again corresponds to a 
spectral flow of surface states similar to class A, but now V ea(fc||) £ a{'H{k\\)), there exists 
^0i~k\\) £ ^(^(“^ 11 )) such that e^(—A:||) = —ea{k\\). Thus for each zero mode at /cy, there 
must also be one at — fey, as shown in subfigure 3II. Hence the parity of the spectral flow in 
class D depends on the existence of an odd number of zero modes at one of the involution 
fixed-points. 

The Z 2 -factor was discussed previously by Kitaev [2] in the context of weak topological 
insulators. We first choose one of the involution fixed-points, then if there is an even/odd 
number of zero modes at the chosen point, we can determine if there is an even/odd number 
of zero modes at the other hxed-point by computing the parity of the spectral flow. In 
subfigure 3III we represent a state with no net spectral flow, that has two zero modes, one 
at each hxed-point. This is the non-trivial 0 01 state in Z 0 Z 2 . Let us remark that there 
is no canonical way of dehning this invariant and a choice of a hxed-point must be taken. 
We also comment that the stability of this class depends on 1) respecting particle-hole 
symmetry and 2) the existence of 2 hxed-points on Xg. In particular 2) does not hold for 
systems where there are umklapp processes [15], leading to phenomena such as a Peierls’s 
transition [48] on the boundary of our half-space. This makes the Z 2 -factor correspond to 
so-called weak topological stability [49]. We briehy discuss the connections of our work to 
weak topological phases in subsection 5.3. 
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TlZ 

Symmetry 

0 


1 

d 

2 

3 

A 

0 

0 

0 

z 

Z 2 

AI 

I 

0 

0 

0 

0 

AH 

-I 

0 

0 

Z2 

Z 3 

D 

0 

I 

Z2 

Z © Z2 

Z^ © Z2 

C 

0 

-I 

0 

z 

Z 2 


Table 1: Classes of topologically stable Fermi surfaces for half-space systems with a 
gapped bulk in dimension d = 1, 2 and 3. 0 is the time-reversal operator and H is the 
particle-hole symmetry operator. We use the names given in [26] (AZ). 0 on the 
right-side denotes the trivial group. 0 on the left-side indicates absence of said symmetry. 
When 0 = H = 0 the surface Brillouin zone is a regular torus and for either 
0^ = ±7 or = ±I it is a torus with involution Higher-dimensional classes relate 

to weak bulk topological phases. 


Examples of systems with a non-trivial class on the Z-factor are spinless p-wave super¬ 
conductors [38] , where Bogoliubov surface zero modes at one of the involution fixed-points 
{—1,1}, satisfy an analogous reality condition to that of Majorana fermions, hence they are 
called Majorana zero modes. Bogoliubov quasi-particles are effectively charge-less, thus, 
our zero modes only contribute to the thermal current. Such a current could be induced 
by a thermal gradient on our system, similar to the electric field for class A [50]. This 
spectral flow has also an anomaly inflow interpretation, where the anomaly is analogous to 
a perturbative gravitational anomaly [50]. 

For symmetry class C, using eq. (4.13) we have 

KH-\S^'^) ^ KO-^ © KO-^ 

RiZ. (5.4) 

Since now = —I, Fermi points at the involution fixed-points must be doubly degenerate 
and the mod 2 invariant must always vanish. This is similar to Kramer’s theorem for time- 
reversal symmetry. Thus the spectral flow must come in pairs and there are no Majorana 
zero modes. There is the same perturbative gravitational anomaly inflow interpretation 
as in symmetry class D (subfigure 3II) for the Bogoliubov zero modes. Condensed matter 
systems with r? = —I in d = 2 correspond to so-called d-wave superconductors [38]. 

For symmetry classes AI and AH the spectral flow must vanish because of the sym¬ 
metry of the spectrum e^(—/cy) = ea{k\\) for some a,j3. By computing class AI., d = 2 

KR-\S^^^) « KO-'^ © KO-^ 

= 0. (5.5) 

We see there is only the trivial class of Fermi points, i.e. all Fermi points are adiabatically 
gappable in class AI, d = 2. For symmetry class AH, d = 2 we obtain 

^ Z2 . (5.6) 
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Due to Kramer’s theorem, if Fermi points arise at the involution fixed-points {1,-1} of X^, 
they must be doubly degenerate, and thus the invariant is the number of double-crossed 
Fermi points mod2 [38]. Such double-crossed Fermi points have a U{l)-spin anomaly inflow 
interpretation [6] and are often called helical Dirac zero modes because of the analogous 
role to the helicity of relativistic Dirac fermions played by spin-orbit coupling. A physical 
example is the quantum spin Hall effect (QSHE) [38]. 

5.3 Higher dimensional cases and weak topological phases 

Our results for all dimensions and respective symmetry classes are summarized in Table 
1. We only point to [43], section 3 and [13], Theorem 11.8 where the methods for com¬ 
puting these groups are presented with the necessary amount of detail. For dimensions 2 
and 3 Table 1 shows non-trivial classes. One could worry that, using the bulk-boundary 
correspondence, these classes should appear in [38] and they do not. We note that we 
have assumed discrete translation-invariance, thus the proper correspondence would be be¬ 
tween Table 1 and weak topological phases [2],[38], that is, those which depend on discrete 
translation-invariance. These are usually constructed by stacking d-dimensional systems 
on top of each other, forming a d-|- 1-dimensional system and looking for gapless edge exci¬ 
tations [38],[49]. We also note that computions of weak topological phases sometimes yield 
more non-trivial classes than those which appear in Table 1. We believe this is due to the 
fact that our half-spaces have less boundaries than the prototypical construction employed 
to compute these, nonetheless a more detailed examination is needed to understand the 
connection between the two. We shall postpone the physical interpretation, possibly of 
anomaly inflow, of the higher dimensional cases for future work. 

6 Conclusions 

We have rigorously derived a classification of topologically stable Fermi surfaces from elec¬ 
tronic band theory and adiabatic evolution for systems on an infinite crystal and systems 
on a half-space. We rigorously proved that there can only be globally topologically unsta¬ 
ble Fermi surfaces for all symmetry classes considered for systems on an infinite crystal. It 
remains to be seen to which, if not all, of the symmetry classes this result extends to. 

We also adapted a simplified version of Hofava’s local stability analysis of Fermi sur¬ 
faces to the single-particle formalism and explained its relation to our global stability anal¬ 
ysis. For systems on a half-space, the construction presented above yields the A-groups 
A-i(X^), Afi-^(Xs), KH-^{Xs),KR-'^{Xs) and KH-\Xs), with X^ the surface Brillouin 
zone, for the cases of discrete translation-invariance and the two different implementations 
of both particle-hole and time-reversal symmetry (classes A, D, C, AI and All). We also 
computed all possible classes for X^ = and d = 1,2,3, summarized in Table 

1. We found that Fermi surfaces belonging to a non-trivial class have an anomaly inflow 
interpretation, with corresponding anomaly charge/spin chiral and gravitational anoma¬ 
lies [6],[7], [50]. The nature of higher dimensional classes presented in Table 1 seems to 
correspond to weak topological phases, but a more detailed analysis is needed. 

We mention in passing that our results for half-space systems strengthen the connection 
found by Hofava [1] between Fermi surfaces and so-called D-branes in Type IIA and Type 
I'A string theory [51], where Hofava obtained the same A-groups but of a compact space- 
time instead of the Brillouin zone. 
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A crucial point clarified in this work is that the allowance of a continuous spectrum for 
the Bloch Hamiltonian 1-L{k\\) by the boundary of the half-space system permits a rigorous 
topological interpretation of non-trivial Fermi surfaces, without resorting to low-energy 
limits. 

Finally, it is necessary to develop a twisted formalism, so that we can include all 
possible symmetries as in [13], [14], overcoming the apparent different choice between the 
norm and the compact-open topology and, if possible, to extend it so that in a single 
mathematical framework, we can encompass all topological phenomena of non-interacting 
fermions in condensed matter. 

A The Riesz metric for unbounded operators and Bloch’s theorem 

We have employed the machinery of [16],[18] to derive a classification of globally stable 
Fermi surfaces and we have imposed that our Bloch Hamiltonian operator %{k) (or 1-L{k\\) 
respectively) be bounded. At first glance this imposition seems unjustified as equations 
(2.2), (2.19) clearly represent unbounded operators. Nonetheless, let denote the 

set of not necessarily unbounded self-adjoint Fredholm operators. We define the Riesz 
metric on CJ^^“(=^) as 

L ■ , (A.l) 

n^nii + n^n)-K (a.2) 

dR : X C.F"“(^) ^ M, (A.3) 

dR{nun2) ^ \\L{ni) - L{n2)l (a.4) 

where jj. jj denotes the standard operator norm and L is known as the Riesz transform [17]. 
The topology on induced by dR is known as the Riesz topology and the natural 

inclusion ^ is a homotopy equivalence. Hence splits into 

3 path-components and is null-homotopic while is again a classifying 

space for the A“^-functor [17]. The issue with the choice of the Riesz topology is that it is 
generally hard to determine whether a family of operators in is continuous [17]. 

However, as is shown in [25] for the infinite crystal, 1-L{k) with H(r) bounded has a fixed 
domain for all G X and as explained in [52], families of operators with a fixed domain 
are continuous in the Riesz topology. Similar arguments apply to the half-space Bloch 
Hamiltonian T-L{k\\) as long as we fix the boundary conditions or, more generally, restrict 
to boundary conditions that have a fixed domain in 
Given that 

L(S'R) = ^LiiU ), (A.5) 

L{Qn) = QL{n) , (A.6) 

the equivariant results considered in section 4 extend to the Riesz topology for CJ^(J^), 
CT{J^q), CTI°‘{J^r) and CFI°‘{J^q) since the corresponding Bloch Hamiltonians and 
boundary conditions have a fixed domain in J^r or respectively. 
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